Helium-4 Synthesis in an Anisotropic Universe 
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Abstract 

We calculate the *He abundance in a universe of Bianchi type 1 whose cosmic anisotropy is dynami- 
cally generated by a fluid with anisotropic equation of state. Requiring that the relative variation of 
mass fraction of *He is less than 4% with respect to the standard isotropic case to be consistent with 
astrophysical data, we constrain the parameter of cosmic anisotropy, the shear S, as |E(rf)| < 0.4, 
where Tf is the freeze-out temperature of the weak interactions that interconvert neutrons and pro- 
tons. Anisotropic fluids, whose energy density is subdominant with respect to the energy content of 
the Universe during inflation and radiation era, generate much smaller shears at the time of freeze- 
out and then do not appreciably affect the standard ''He production. This is the case of anisotropic 
dark energy, and of a uniform magnetic field with energy density much smaller than about 1.25 
times the energy density of neutrinos. 

PACS numbers: 98.80.Jk, 98.80.Ft 



I. INTRODUCTION 

The high level of isotropy of the cosmic microwave 
background (CMB) radiation is the most convincing 
justification of the Cosmological Principle: the Uni- 
verse is homogeneous and isotropic at large cosmological 
scales However, tiny deviations from perfect isotropy 
are not excluded by present CMB data. Indeed, a par- 
ticular anisotropic cosmological model of Bianchi type 
I, known as ellipsoidal universe [2I, Isj, can better match 
CMB data and solve the so-called "quadrupole problem," 
namely the lack of CMB power detected on large angular 
scales. 

Various mechanisms could give rise to an ellipsoidal 
universe, such as a uniform cosmological magnetic 
field [3tJ1j topological defects (e.g. cosmic stings, domain 
walls) [4|, or a dark energy fluid with anisotropic equa- 
tion of state 0, 0| • Independently on the nature of the 
mechanism, however, a modification of the standard pic- 
ture of primordial nucleosynthesis can occur if universe 
anisotropization takes place during the early Universe 

The aim of this paper is, indeed, to constrain the level 
of cosmic anisotropy, so as to be consistent with obser- 
vational bounds on primeval "^He abundance. 



II. ELLIPSOIDAL UNIVERSE 



The ellipsoidal universe is a cosmological 

model described by the Bianchi I line element [l^ 



ds^ = df - a^{t){dx^ + dy^) - b^{t) dz^ 



(1) 



with two scale factors, a and 6, normalized as a = 6 = 1 at 
the present cosmic time. Cosmic anisotropy is quantified 



by the shear 



Ha — Hb 
2Ha + Hb 



(2) 



with Ha = a/a and Hb = b/b, while H = A/A = {2Ha + 
Hb)/3 and A = {a'^by/^ play the usual role of Hubble 
and expansion parameters, respectively. (Here and in 
the following a dot indicates a derivative with respect to 
cosmic time t). 

Anisotropy of the Universe is not assumed a priori but 
dynamically generated by an anisotropic fluid (A) with 



two equations of state: w\ = p\/pA and m\ = p^/pA, 
where and are respectively the pressures along the 
X (y) and z directions, and pA the energy density. The 
source of cosmic anisotropy is then parameterized by the 
skewness Sa = w\ — i(r^, while wa = i2w\ + Wj^)/3 
represents the usual equation of state parameter. 

Friedmann equation in ellipsoidal universe takes the 
form [H [ill 



(1 - j:^)h^ 



SttG 



(p + Pa), 



(3) 



where p is the sum of the energy densities of the usual 
components in the standard model, namely photons Pj, 
neutrinos p^,, matter pm, and dark energy pde- (In the 
following discussion, we neglect the effects of matter since 
nucleosynthesis takes place in radiation dominated era.) 

The shear is sourced by the skewness according to the 
equation [H, |li| 



(i/S)- +3H^^ 



BttG 

IT 



{pijSty + paSa), 



(4) 



where Si, is the neutrino skewness that takes care of ef- 
fects of anisotropic distribution of neutrinos. It depends 
on the shear and its form will be discussed later. 
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Inflation generally causes an isotropization of the Uni- 
verse: any cosmic shear present before and/or during 
inflation will be reduced to a vanishingly small value af- 
ter inflation (see discussion below). Nevertheless, if a 
source of anisotropy is present after inflation (e.g. an 
anisotropic fluid), then the cosmic shear can grow and 
be different from zero at the time of decoupling. If this is 
the case, planar cosmic symmetry induces a quadrupole 
term in the CMB radiation which adds to that caused 
by the inflation-produced gravitational potential at the 
last scattering surface. If the planar-metric induced 
quadrupole is comparable to the inflation-produced one, 
the overall quadrupole power may match the "anoma- 
lously low" value of the observed quadrupole [1, 0] ■ The 
capability to solve the CMB quadrupole problem is the 
main attractive feature of the ellipsoidal universe model. 



III. HELIUM-4 SYNTHESIS 



The mass fraction of ^He produced by standard pri- 



mordial nucleosynthesis at the cosmic time t 



(0) 



300 s 



corresponding to a temperature of 



O.OTMeV 



IS 



lit 



y(0) 



2{n/p) 



(0) 
nuc 



where 



(n/p) 



(0) , 

nuc 



1 + (n/p) 



-Q/Tf 



(0) 

nuc 



0.25, 



1/7 



(5) 



(6) 



is the neutron-to-proton number density ratio. [We in- 
dicate quantities in the standard isotropic cosmological 
model with an index "(0)".] The first exponential factor 
in Eq. ([6|), with Q ~ 1.3MeV being the neutron-proton 
mass difference, is the neutron-proton number density ra- 
tio at the time of freeze-out, namely when the expansion 
rate of the Universe, equals the rate for the weak 

interactions, T ~ G^T^, that interconvert neutrons and 
protons (Gf is the Fermi constant). This happens at a 
temperature of about ~ O.SMeV 
terium bottleneck" llj 



17|. Due to "deu- 
the production of ''He is delayed 
until the Universe has cooled to the temperature Tnuc- 
In this time lag, neutrons decay reducing their relative 
abundance and, in turn, that of *He. This gives the sec- 
ond exponential term in Eq. (jS]), where r„ — 885 s is the 
mean neutron lifetime. 

In ellipsoidal universe, both the freeze-out temperature 
and the time of nucleosynthesis are modified, and so is 
''He abundance. 

Astrophysical observations fix the value of ^He abun- 
dance as ~ 0.25 ±0.01 [i3j. (See ^ and references 
therein for more recent estimates of Y^^^^ which are, how- 
ever, all consistent with that quoted here. This can be 
considered as the most conservative estimate of since 



it possesses the largest uncertainty.) Therefore, to be 
consistent with experimental data, we must require that 
the maximum variation of ^He abundance (with respect 
to the isotropic case) is below the 4%. 

A general expression for the freeze-out temperature in 
ellipsoidal universe is easily obtained from Eq. ([3]) if one 
assumes that the energy content of the anisotropic fluid 
is subdominant with respect to that of radiation: 



Tf 



T 



(0) 



/ 



(1 



I]2)l/6 ' 



(7) 



where S/ is the shear at the time of freeze-out and we 
used the fact that p = {n^ / 30) g^,T'^, with g^, the total 
number of effectively massless degree of freedom [l3| • In 
the following we simply assume 5, = 3.36 during nucle- 
osynthesis (even if g^, increases from that value to 10.75 
near Tf). 

The time when ''He is produced is easily found by in- 
tegrating the Friedmann equation with respect to time: 



3V5 



27r3/2 



1/2 



-(1 



S^) 



2n1/2 



(8) 



where mpi is the Planck mass and we used the fact that 

Due to positivity of the energy and looking at the 
Friedmann equation we see that the shear is bounded 
in the interval [—1,11. This implies, using Eqs. ^ and 



that Tf > rl°^ and t 



_^ ''nuc- 



Since the He abun- 



with rj"'' and 



dance, Y, is given by Eqs. ([5|) and 

tnuc replaced by Tf and tnuc, we conclude that in ellip- 
soidal universe there is an overproduction of *He with 
respect to the isotropic case, whatever is the nature of 
the anisotropic source. 

In order to calculate this positive variation of ^He mass 
fraction, one has to specify the anisotropic source so as 
to integrate Eq. ([4]), flnd the shear as a function of tem- 
perature, and obtain tnuc from Eq. ([5]). 

This can be done analytically only in the case where 
the effects of the skewness are neglected {6a — Si, = 0). 
Indeed, the case = is that studied numerically in the 
literature taking into account both the full set of nuclear 
reactions leading to the production of light elements and 
the effects of anisotropic distribution of neutrinos 2^, 21 1 . 
The effects of having 5^ 7^ are studied below, but we 
will show that they are negligible (at least for small values 
of the shear) . 

Introducing the anisotropy parameter B — ^^/(l — S^), 
the shear equation ^ gives B (x T^, so we can easily 
solve Eq. ([7]) for the freeze-out temperature, and perform 
the integral in Eq. ([5]) to get the time of nucleosynthesis: 



T;=T)"Vi[S(r)"0], 

tnuc = tiuc f2[B{Tnuc)] 



(9) 
(10) 
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FIG. 1: Relative increase of He abundance in ellipsoidal 
universe with zero skewness Sa and neglecting neutrino 
anisotropy effects (continuous line) as a function of the 
anisotropy parameter Bo = B{T = 50 x lO^K), where 
B — — Tt^). The dotted line is the asymptotic expan- 

sion Bo X 1%, while 300% is the limiting value for Bo ^ oo. 
The dashed line is the relative increase of *He abundance in 
the same cosmological model but assuming no variation on 



the time of nucleosynthesis, tn 



where 



(0) 




2 X 31/3 + 21/3(9 + V81 - 122)3)2/3 



62/3(9 + V81- 12x3)1/3 
— X arccoshv^ . 



(11) 
(12) 

In Fig. 1, we plot the relative increase of ^He abundance 
(with respect to the isotropic case) as a function of Bq, 
namely the anisotropy parameter evaluated at the refer- 
ence temperature Tq = 50 x lO^K, well before the nucle- 
osynthesis starts. (B, and then Bq, are the same quantity 
introduced in [2l|.) The asymptotic expansions of such 
an increase, for small and large values of Bq, are: 



y(o) 



Bo 
Bo 



0, 

00, 



(13) 



where 



y(0))Qy(0) 



12 r2 



1 - y(°) 
y(o) 



C'(lnBo) ~ 0.01, (14) 



(15) 



The first term in the right hand side of Eq. takes 
into account the rise of the freeze-out temperature in 
anisotropic universe, while the logarithmic term takes 
care of the reduction of time lag between the freeze-out 
and the end of nucleosynthesis, and is negligible with re- 
spect to the first one. 



The numerical analysis of [21| shows an increase of 
light element abundances. In particular, the relative in- 
crease found for ^He is linear for moderate values of Bn 



{Bq < 10) and is about Bo x 3%. Therefore, our over- 
simplified analytical model confirms qualitatively (and 
to some extent also quantitatively) the numerical results 
of 21]. 1 

The dashed line in Fig. 1 is the relative increase of ''He 
abundance assuming no variation on the time of nucle- 
osynthesis, tnuc = inuc- As it is clcar from the figure, the 
time delay efi'ect due to cosmic anisotropy, tauc ^ inuc, 
causes appreciable effects only for large shears which are, 
however, unrealistic because of the large excess of '^He. 
We note that if we just replace Tj"' with Tf and leave 



in Eqs. dSD and (ED, we obtain a lower limit on Y. 
Imposing that the mass fraction of "iHe is less than 4% 
with respect to the standard isotropic case, we obtain a 
conservative, but model-independent limit (not depend- 
ing on 6a and wa) on the level of cosmic anisotropy at 
the time of freeze-out: 



|S(T/)|<0.4. 



(16) 



It is straightforward to show that the above limit is in 
agreement with the limit obtained by translating the cur- 
rent bound on the total number of effectively massless 
degree of freedom at time of freeze-out. Indeed, assum- 
ing as before that the energy content of the anisotropic 
fiuid is subdominant with respect to that of radiation, 
we can rewrite Eq. ([3]) as the usual Friedmann equation 
= (87rG'/3)/9 where now p = (7r2/30)g,,cff T"* with 



J*,cS — 



1-E2 



(17) 



Therefore, the effect of having a nonzero shear at the 
time of freeze-out can be regarded as a change in the 
total number of effectively massless degree of freedom, 
which is usually parameterized by the effective number 
of neutrino species, N^, as 17 1 



11 7^/4 



4/3 



(18) 



(The standard value ~ 10.75 near Tj corresponds to 
take = 3 in the above equation.) Using the cur- 
rent bound = 3.2 ± 1.2 (95% C.L.) [19] on the ef- 
fective number of neutrino species at the time of freeze- 
out, and comparing Eqs. PT)) and we get |I](T/)| = 
0.11 ± 0.34 (95% C.L.), where we used the Gauss error 
propagation law to propagate the uncertainty on Ni, to 



1 It is worth noticing that our analysis needs to be modified for 
very large values of the anisotropy parameter B sinee, as pointed 
out in [23 , the equilibrium of weak interactions can be broken by 
very high levels of anisotropic expansion. However, our forthcom- 
ing results will rely just on the part of Fig. 1 that corresponds to 
moderate values of the anisotropy parameter (namely Bo S, 10). 
whose validity has been already confirmed numerically in [2l| . 
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|S(T/)|. So, we obtain the upper bound |E(r/)| < 0.45, 
which is compatible with Eq. (jl6p . 

The effects of anisotropic distribution of neutrinos can 
be described as follows. For temperature greater than 
t'^^^ = ©(MeV), neutrinos are strongly coupled to pri- 
mordial plasma, so their distribution is isotropic and no 
neutrino skewness results. Below a temperature slightly 
lower than Tj'^\ instead, neutrinos begin to free-stream 



and generate a skewness 



dT' 



(19) 



Here, for the sake of simplicity, we assumed an instan- 
taneous neutrino decoupling at rj'^^ ~ Tf, so neutrino 
free-streaming affects only the time when "^He is pro- 
duced, namely t^uc- The above result is valid for 
small values of the shear (|E| <C 1) and can be obtained 
from [20| taking the limit, in the Misner's anisotropy po- 
tential, of large collision time {tc — ^ oo) for the typical 
reactions of neutrinos with plasma. 

Using (HH) in Eq. (g]), we get for T <Tf 



1/2 



cos 



,ln^' 

Tf. 



1 . 

— sm 



Cu In— 

(20) 

where = ^8n,^/5 — 1/4 with fl^ — p^/pa being the 
neutrino energy density parameter and pa — 3H^/8ttG 
the critical density. The neutrino energy density param- 
eter riy is constant during radiation era and, assuming 
three neutrino massless species, equal to about 0.4 after 
neutrino decoupling [l7| . so Ci, ~ 0.6. The shear is an 
oscillating function of time with a damping factor pro- 
portional to T^^/^. This leads to very tiny variations of 
^nuc with respect to the case 6^ = 0, and gives rise to a 
small increase of {Y -Y^°^)/Y'^°'> x 100%, which is below 
the 0.15% for < Bo < 1 (corresponding to S < 0.3). 

Since the absolute value of the shear at the time of 
freeze-out must be significatively smaller than one [see 
Eq. pop ] , we can now consider a simplified (but more re- 
alistic) model where the shear is a small quantity during 
radiation era and the effects of both neutrino skewness 
6i, and external anisotropic sources 6a are taken into ac- 
count. 

For small shears and 6a constant, the energy density 
of an anisotropic fluid evolves as in the case of isotropic 
universe, pA oc ^-3(1+"' ^i) 14 1^ and the shear equa- 
tion (|4]) can be solved to give in radiation era and before 
neutrino decoupling. 



constant 
A 



A,0 



2 — 3wa ^r.O 



A 



1 — 3WA 



(21) 



where we assumed that the energy density of the 
anisotropic component is small with respect to that of 
radiation. This is the same as assuming wa < 1/3, or 



6a^a,o ^ ^r,o if WA = 1/3. Here, ilA,o and Qr,o are 
the present energy density parameters of anisotropic fluid 
and radiation, respectively. ^ 

We can fix the integration constant in Eq. ([21]) by 
evaluating the shear at early times, for example at the 
end of inflation, A = A^nd ^1. If wa < 1/3 we get 
constant ~ ^end^^end, whcrc Send = 2](Aond), whilc for 
WA = 1/3 we have constant ~ Acndi^cnd — 6A^A,o/^r,o)- 
As shown in 0] , any cosmic anisotropy is washed out (ex- 
ponentially) during (de Sitter) inflation (for wa > — 1), 
so that anisotropy can develop just at the end of inflation 
starting from a vanishingly small value. ^ This means 
that Send — 0, so that we can neglect the first term in 
the right hand side of Eq. ([2T|) for A > ^end- 

Let us assume for the moment that the effects of neu- 
trino free-streaming are negligible. In this case, the above 
solution ((2T|) is valid throughout nucleosynthesis and, 
since E <C 1, we conclude that no appreciable changes 
on '^He production occur with respect to the isotropic 
case. We can now check the validity of the assumption of 
neglecting neutrino anisotropy. By inserting Eq. (1211) in 
Eq. (fT9)) we find that the ratio of the anisotropy sources 
in Eq. (g]) is, for T <Tf and wa ^ 1/3: 



Pi^6u 
Pa6a 



8 r!,,o 1 - [T/TfY- 



5 VLrfl (1 — 3wa)(2 — 3wa) 



(22) 



where J7y_o is the present neutrino energy density pa- 
rameter. Assuming three neutrino massless species, we 
have rijy.o/f^r.o — 0.4 17|. For the cosmologically inter- 
esting cases of anisotropic dark energy {wa — ~l)j a 
cosmic domain wall {wa — ~2/3), and a cosmic string 
{wa = ^1/3), the absolute value of the ratio ([22]) is max- 
imum at r = T„uc and is much smaller than one (0.03, 
0.05, and 0.11, respectively, assuming Tf ~ T'j^''), and 
this justifies our previous assumption. 

The case wa — 1/3, namely an anisotropic compo- 
nent of radiation type, has to be analyzed separately. 
In this case, the shear equation (|4]) can be solved and 
gives, for T < Tf, Eq. ^ with the factor 1/c^ multi- 
plying the sine function replaced by ca/c^, where ca = 



^ For Wa = 1/3, Eq. II21I I is correct up to a logarithmic 
term. Indeed, as shown in Q], the last term in the right 
hand side of Eq. I I21I I should be divided, in this case, by 
1 + 25^ (Oyl_o/f^r,o) lri(^/^ond)i where A^nd is the expansion pa- 
rameter at the end of inflation. However, the inclusion of this 
term modifies Eq. I I21I I only to second order in the small quantity 
SA^A,o/^r,o- Therefore, for simplicity, we neglect this term in 
the following. 

In de Sitter inflation, subdominant anisotropic fluids are such 
that Wa > —1, or 5^f24 <S 1 if wa = —1- For such fluids 
and in the limit of small shears, it is easy to see that Send — 
Sie-3iv if > and E^nd ^ -((5a/3«)a) 0^,0 e-3(i+™>i)JV 
if — 1 < WA < and A 2> A^. Here, is the shear at the 
beginning of inflation at A = Ai, and A'^ > 60 the number of 
e- folds of inflation since inflation began [l7tl . 
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Sa^a/^! - 1/2. Using Eq. (HH) evaluated at T = T/, 
we find ca = ^r,o^A/^A,o — 1/2. Since both ^a and 
fir are constant in radiation era and scale as T in mat- 
ter era, we get ca — il^ ~ 1/2 in radiation era, where 
— ^ — — 0.6 after neutrino decoupling. Therefore 
CA — 0.1. Since neglecting neutrino anisotropy we found 
that S is constant (up to a logarithmic correction) in 
radiation era [see Eq. ([2T|) ] and does not affect Helium-4 
synthesis, we conclude that also in the case (5y 7^ 0, where 
the shear is an oscillating function of time with a damp- 
ing factor proportional to T"^/^, no appreciable changes 
on ^He production occur with respect to the isotropic 
case. 

Before concluding, let us include in our analysis a com- 
ponent of free-streaming gravitons, for inflation generally 
predicts gravitational waves, namely tensor fluctuations, 
which are not in thermal equilibrium below the Planck 
scale and then can be considered as free-streaming radia- 
tion from inflation until today. Gravitational waves intro- 
duce the extra term (87rG/3)pGW<5GW on the right-hand- 
side of Eq. dU, where pcw and (5gw are the graviton en- 
ergy density and skewness, respectively. The energy den- 
sity associated to this background of gravitational waves 
is typically very small: i^Gw.o/^^r.o ^ 10^^ for modes 
that cross inside the horizon while the Universe is radia- 



tion dominated [17|, where ricw.o is the present energy 
density parameter of gravitons. The graviton skewness is 
given, after inflation, by an expression similar to Eq. (jl9p 



3gw 



(T) 



S(T'), 



(23) 



where Trh is the so-called reheating temperature, that is 
the temperature of the cosmic plasma at the beginning of 
the radiation era. (Here and in the following we assume 
that the reheating phase, during which the energy of the 
inflaton is converted into ordinary matter is "instanta- 
neous" so that, after inflation, the universe enters directly 
into the radiation era.) Now it is easy to show that, 
due to the smallness of f^GW,o, the effect of gravitational 
waves in the evolution of the shear is completely negligi- 
ble. In fact, proceeding as we did in obtaining Eq. (1^ . 
we can verify that the ratio \pgw5gw I Pa5a\ is much 
smaller than unity. Indeed, for T < Trh and wa ^1/3, 
it is given by the right hand side of Eq. (|22|) with V,^^ 
and Tf replaced by f^ow.o and Trh, respectively. There- 
fore we have \pgwSgw/paSa\ ~ ^GW,o/^r,o ^ 1. For 
T < Trh and wa = 1/3 we get, instead, 

Pgw^gw _ 8 r^Gw.o 
^5 ~ 



ln(TRH/T). 



(24) 



paSa 5 rirfi 

The absolute value of the above ratio is maximum for 
T = Tnuc and for the largest allowed value of Trh, Trh — 
10^''GeV [13]. Also in this case it is much smaller than 
unity: \pgwSgw / Pa5a\ ^ TSilcwfi/^rfi < 1. 

Let us conclude by observing that for a uniform mag- 
netic field WA = 1/3 and 6a = 2. Therefore, the above 



results show that uniform magnetic fields created at in- 
flation and whose energy density is small with respect to 
that of radiation do not affect nucleosynthesis. However, 
in the presence of an external uniform magnetic field, 
nucleosynthesis is affected, other than by the effect of 
anisotropization of the Universe due to a nonvanishing 
shear, also by the increase of weak reaction rates, of the 
expansion rate of the Universe, and of the electron den- 
sity [23|. Taking into account all these effects, but not 
the effect here studied of nonvanishing E, the authors 
of 23] found that observations of light elements are com- 
patible with a magnetic field energy density lower than 
PB ^ 0.28p^, where pB = is the magnetic energy 

density associated to a uniform magnetic field of intensity 
B. Their analysis is correct as long as the effect of the 
shear can be neglected which means, in light of the previ- 
ous discussion, that the magnetic field must be a subdom- 
inant component of the Universe during nucleosynthesis. 
This is indeed the case, since the subdominance condition 
for a uniform magnetic field, 2r2g.o ^ ^r,0: translates to 



„ ^ max 

Pb^ Pb - 



1.25p^ 



(25) 



after neutrino decoupling, a limit about 5 times greater 
than that allowed by the analysis of _23] . 

It is worth noticing that the above limit on the inten- 
sity of a cosmological magnetic field is much less stringent 
than that coming from the analysis of the CMB radia- 



tion @,[l[ii,[2i 



25| which is at least two order of magni- 
tude stronger. This agrees with Barrow's result Q that 
anisotropic fluids that create temperature anisotropics 
compatible with CMB spectrum do not have a signifi- 
cant effect on the primordial synthesis of ''He. 



IV. CONCLUSIONS 

In this paper, we have analyzed the effects caused by 
cosmic anisotropy on the primordial production of '^He. 
We worked in the context of a cosmological model of 
Bianchi type I where the anisotropy of spatial geome- 
try, the shear E, is generated by a fluid with anisotropic 
equation of state. 

We found that in such an anisotropic universe there 
is an overproduction of ''He with respect to the stan- 
dard isotropic case. Imposing that the relative increase 
of ^He abundance is below the 4% to be consistent with 
observational data, we constrained the absolute value of 
the shear to be less than 0.4 at the time of freeze-out. 
This limit does not depend on the equation(s) of state 
of the anisotropic fluid and has been obtained assuming 
that the energy density of the anisotropic fluid is small 
compared to that of radiation. 

Moreover, we showed that anisotropic fluids gener- 
ated at inflation, such as dark energy with anisotropic 
equation of state and a uniform magnetic field, create 
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anisotropies much smaller than the above limit if their 
energy densities are subdominant with respect to that of 
the Universe during inflation and radiation era. In par- 
ticular, the existence of a uniform magnetic field at the 
time of nucleosynthesis is compatible with astrophysical 
data if its energy density is much smaller than about 1.25 
times the energy density of neutrinos. 

We would like to thank P. Cea and J. D. Barrow for 
very helpful discussions. 
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